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Abstract. The scattering matrix of the Schrodinger operator with smooth short-range 
electric and magnetic potentials is considered. The asymptotic density of the eigenvalues of 
this scattering matrix in the high energy regime is determined. An explicit formula for this 
density is given. This formula involves only the magnetic vector-potential. 



1. Main result and discussion 

1.1. Main result. Consider the Schrodinger operator 

H=(iV + A) 2 + V in L 2 (R d ), d>2; 

here V : M. d — > R is an electric potential and A = (A%, . . . , Ad) : M d — > M. d is a magnetic 
vector-potential. We assume that both V and A are infinitely differentiable and satisfy the 
estimates 

(1.1) \d a V(x)\ < C a {x)-"-W \d a A(x)\ < C a {x)-"-W, p>l, 

for all multi-indices a; here (x) = (1 + |x| 2 ) 1//2 . Let H = —A; we denote by S(k) the 
scattering matrix associated with the pair H, H at the energy A = k 2 > 0. We recall the 
precise definition of the scattering matrix in Section I2.1( here we only note that S(k) is a 
unitary operator in L 2 (§ d_1 ) and the operator S(k) — I is compact (see e.g. [TH HE]). Thus, 
the spectrum of S(k) consists of eigenvalues on the unit circle T; all eigenvalues (apart from 
possibly 1) have finite multiplicities and can accumulate only to 1. Our aim is to describe 
the asymptotic behaviour of these eigenvalues as k — > oo. 
If A = 0, one has 

(1.2) \\S(k) =0{k~ 1 ), k->oo (A = 0), 

and so the eigenvalues of S(k) are located on an arc near 1 of length O^ -1 ). The large energy 
asymptotics of S(k) in this case is given by the Born approximation, see e.g. [TH Chapter 8]; 
this makes the analysis of S(k) rather explicit. In [6], using the Born approximation, we 
have determined the large energy asymptotic density of the spectrum of S(k) for A = 0; we 
will say more about this in the next subsection. 

When A ^ 0, the situation is radically different: as k — > oo, the norm ||if>(A;) — J|| does not 
tend to zero and the Born approximation is no longer valid. The high energy asymptotic 
expansion of the scattering amplitude (= the integral kernel of S(k) — I) for the magnetic 
Schrodinger operator was obtained (through a very difficult proof) in [17] . Our main result 
below is merely a spectral consequence of |17j . 
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We need some notation. For any u G S d_1 , let C M. d be the hyperplane passing through 
the origin and orthogonal to u>: 

A LU = {x eR d : (x,u) = 0}. 

We equip both §> d ~ l and A w with the standard (d — l)-dimensional Lebesgue measure. We 
set 

/DC 
(A(tu + £),u)dt, we§" (eR d 
-oo 

Our main result is as follows: 

Theorem 1.1. Let V , A satisfy (11. ip . Then for any function if G C(T) that vanishes in a 
neighbourhood of the point 1, 

(1.4) Mm k- d+1 Trip(S(k)) = (2n)- d+1 [ [ ip(e iM{uj ^)d£duj. 

Under the assumptions of Theorem ll.il the operator ip(S(k)) has a finite rank, and so the 
trace in the l.h.s. exists; also, the integrand in the r.h.s. vanishes for all sufficiently large 
|£|, and so the integral is absolutely convergent. 

The striking feature of formula (jl.4p is that its r.h.s. does not depend on the electric 
potential V. 

1.2. Discussion. Weak convergence of measures: Theorem 11.11 can be rephrased in terms 
of weak convergence of the eigenvalue counting measures. Let us denote the eigenvalues of 
S(k) by e l9n ( k \ where n G N and 9 n (k) G [— tt, tt); as usual, the eigenvalues are counted with 
multiplicities taken into account. Let i C T \ {1} be any open arc separated away from 1; 
we set 

(1.5) /i k (t) = #{n G N : e ie " {k) G t}, k > 0, 

(1.6) MO = (2vr)- d+1 / / d£du,, IL(0 = {( e A u : e iM ^ G l}, 

where # represents the number of elements in a set. The measures fi can be alternatively 
defined by requiring that 



;i.7) Tr<p(S(k)) = / <p{z)dii k {z) 

JT 

;i.8) (2n)- d+1 [ [ ip(e m ^)d^du = [ f(z)d^z), 



for any continuous function ip on T vanishing in a neighbourhood of 1. With this notation, 
Theorem 11.11 can be rephrased as the weak convergence 

k~ d+1 fik ->■ — y oo. 
The measure /z may be singular at 1, i.e. /i(T \ {1}) may be infinite, but 

(1.9) I \z- l\ e dn{z) < +oo, W > (d- l)/(p- 1). 



The relation (II .9p follows from the elementary estimate on M(u;,£) which is a direct conse- 
quence of (II. ip : 

(1.10) |A^(o;,e)l < C(A)(1 + ICI) 1 ^, Loe§ d -\ £GA W . 
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The case A = 0: In this case, the limiting measure /i vanishes, and in order to obtain 
non-trivial asymptotics of /i^, the problem requires an appropriate rescaling. By (11. 2p . the 
spectrum of S(k) consists of eigenvalues which lie on an arc of length 0(fc _1 ) around 1. This 
suggests the following rescaled version of the problem: for an interval S C R \ {0} separated 
away from zero, set 

= #{n G N : kO n {k) e 5}. 

Then it turns out (see [Bj) that 

(1.11) k~ d+1 Jl k — > Jl weakly as k — > oo, 

where the measure Jl is defined as follows. Let 

1 f°° 

(1.12) X(u},£) = -- V{tuj + Z)dt, tueS d -\ eeA, ; 



then 

Jl(5) = (27r)- d+1 [ [ d£du, ft (i) = {( 6 A u : X(u, f) G 5}, 
where 5 C R\{0}. 

Two and three dimensional cases: Let us discuss assumptions (II. ip in dimensions d = 2,3. 
In dimension d = 3, a magnetic vector-potential A satisfying ( II. ip can be constructed for 
any smooth magnetic field B : R 3 — > R 3 (B = curl A) such that div B = and 

(1.13) \d a B(x)\ < C a (x)- p - x -\ a \, p > 1. 

In dimension d = 2, the magnetic field 5 : R 2 — > R, in addition to ( 11. 13ft . must satisfy the 
zero flux condition 

(1.14) $ = / fi(x)dx = 0, B(x) = ^ - 

See [18] for the details of this construction. The zero flux condition is unavoidable in the 
following sense. Let B : R 2 — > R be given which satisfies the estimates (11.131) but the flux 
$ 7^ 0. Then any magnetic vector-potential A(x) for this field will necessarily fail to be 
short-range (i.e. (II. ip fails) but one can construct A(x) with the behaviour ~ 
as |x| —> oo. In this case the scattering theory for H Q , H can still be constructed, but the 
difference S{k) — I will not be compact, see [18] for a detailed discussion and a description of 
the essential spectrum of S(k). A particularly well known example of this is the Aharonov- 
Bohm effect [2]. Thus, in this case the measure fi^ cannot even be defined and the question 
of the spectral asymptotics of the scattering matrix cannot be approached in the same way 
as in Theorem 11.11 

Gauge invariance: The scattering matrix is gauge invariant in the class of short-range 
magnetic vector-potentials. More precisely, if 

(1.15) A(x) = A(x) + V0(x), where (f>(x) — > as |x| — > oo, 

then the scattering matrix S(k) associated with the pair H = (iV + A) 2 + V, H = —A, 
coincides with S(k). See [TTJ [18] for further details and for an interesting discussion of 
examples when S(k) is not gauge invariant (when <p{x) is homogeneous in \x\ of order zero). 
An inspection shows that M(u, £) is gauge invariant under the gauge transformations of the 
class (051) . 

Related work: Although the study of various aspects of asymptotic distribution of eigen- 
values of differential operators has become a well-developed industry, very little attention in 
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mathematical literature has been devoted to analogous questions for the scattering matrix. 
There has been, of course, much work on the "average" characteristics of the scattering 
matrix S(k): on the scattering phase argdet S(k) and on the total scattering cross section; 
but the study of the distribution of individual eigenvalues of S(k) has been less popular. We 
are only aware of mathematical works [SI IH El ED] on this subject. In [5J, the asymptotic 
behaviour of the phases 9 n (k) was determined for a fixed k and n — > oo for potentials V 
with compact support and A = 0. In [I], the same problem was considered in the case of 
potentials V with a power asymptotics at infinity. In [5], this question was studied in the 
presence of a periodic background potential. In [20], the spacing between the phases 9 n {k) 
was studied in a rather special model of scattering on manifolds. 

In the physics community, the question of asymptotic distribution of eigenvalues of the 
scattering matrix has been addressed; see e.g. the works by U. Smilansky and his collabora- 
tors [8] on the eigenvalue statistics of S(k) for obstacle scattering. 

The limiting measure \x arises via integration over straight lines, i.e. over the trajectories 
of the free dynamics, see (11 .3D . Similar in nature asymptotic formulas are known in other 
problems, involving the discrete spectrum of differential and pseudodifferential operators; 



1.3. Key steps of the proof of Theorem ll.il Our proof is heavily based on the results of 
[T7] . In [IT], D. Yafaev suggested a high energy asymptotic expansion for the integral kernel 
of the scattering matrix S(k). This expansion is constructed via approximate scattering 
solutions to the stationary Schrodinger equation 



We recall this construction in Section [2j Using Yafaev's expansion, it is easy to represent 
the scattering matrix S(k) (modulo a negligible error) as a semiclassical pseudodifferential 
operator (\l/DO) on the sphere with the semiclassical parameter h = k~ x and the principal 
symbol e lM ^'^. This representation is already present, in a somewhat different form, in 
[PT] ; we re-derive it in Section [2] in the form convenient for our purposes (see Lemma [2.5p . 
After this, using the standard pseudodifferential techniques, we prove (see Lemma 13.11) the 
asymptotic formula 



for any integers l\ > 0, £2 > such that the sum i\ + £ 2 is sufficiently large. From here it is 
easy to derive the main result by an application of the Weierstrass approximation theorem; 
this is done in Section [3J 

The difference between the case A = and the general case can be understood as follows. 
As mentioned above, the scattering matrix can be approximated by a \[>DO on the sphere 
with a symbol which depends on k. The leading term of the asymptotic expansion of this 
symbol in powers of k~ x is e lM ^ u ^'\ this term involves only the magnetic vector potential A. 
The electric potential turns up only in the next term of the expansion, which (if A = 0) has 
the form ik _1 X(u, £) (see (I1.12p ). Of course, this is related to the fact that the magnetic 
potential is a perturbation of the order 1 (as a differential operator) of —A, whereas the 
electric field is a perturbation of the order 0. 



see pain ma nm. 



Hu = k 2 u. 



(1.16) 
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2. The scattering matrix as a vpDO on the sphere. 

2.1. The scattering matrix. Let us briefly recall the definition of the scattering matrix; 
for the details, we refer to any textbook on scattering theory, e.g. [19] - It is one of the 
fundamental facts of scattering theory that under the conditions fll.ip . the wave operators 

W± = s-lim e itH e- itHo 

t— >±oo 

exist and are complete. The scattering operator S = W+W_ is unitary in L 2 (M. d ) and 
commutes with H . Let F : L 2 (R d ) — > L 2 ((0, oo); L 2 (§ d_1 )) be the unitary operator 

(Fu)(k,u) = -^k( d - 2)/2 u(kcu), k>0, L0EE> d -\ 
V2 

where u is the usual (unitary) Fourier transform of u. The operator F diagonalises Hq, i.e. 

(FH u)(k,uj) = k 2 (Fu)(k,u), Vw G C™(R d ). 

The commutation relation SHq = HqS implies that F also diagonalises S, i.e. 

(FSu)(k, •) — S(k)(Fu)(k, •), 

where S(k) : L 2 (§ d_1 ) — > L 2 (S d_1 ) is the unitary operator known as the (on-shell) scattering 
matrix. 

2.2. Pseudo-differential operators on the sphere. For every u G S d_1 , we identify the 
cotangent space T*E> d ~ l with the plane A u = {i e l d : (^w) = 0} in a standard way. For 
a symbol a G C^°(T*S ci_1 ) and a semiclassical parameter h G (0, 1), the semiclassical ^DO 
Op ft [o~] in L 2 (S d_1 ) is defined via its integral kernel 

(2.1) Op h [a){u,u') = (2nh)- d+l [ e^^'^/V^, £)d£, 

where u, u' G § d_1 . This definition can be extended in a standard way to symbols o satisfying 

(2.2) \d^a(u,0\<C a ?(0~ m ~ lal , weS", £gA„, 

for some m G R and all multi-indices a, We will only be interested in the case m > 0; 
then by the Calderon-Villancourt theorem (see e.g. [13]) combined with a scaling argument, 
Op/Jcr] is bounded and 

(2.3) sup ||Op fc [<7]|| <C(a). 

0<h<l 

Next, if a satisfies (12.21) with m > d — 1, then (see e.g. [121 E]) Op/Jer] is trace class and its 
trace can be computed by integrating the kernel (12.11) over the diagonal: 

(2.4) Tr(Op,[a]) = (2^)- d+1 / / (r(u,£)d£du. 

We will also be interested in symbols which depend on h. For m > 0, let <S((£)~ m ) De the 
class of C°°-smooth symbols a = cr(u>, £,h), h G (0, 1), satisfying the estimate (12. 2p uniformly 
in h G (0, 1) for all multi-indices a, (3. We will need a standard statement about the leading 
term spectral asymptotics of a semiclassical ^DO: 
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Proposition 2.1. Let m > and let o G 5((£)~ m ) be a symbol which admits the represen- 
tation 

(2.5) a = a + ho~\ 

with ao,o~i G 5((£) _m ) and <To independent of h. Then for any non-negative integers £\, £2 
such that £1 + £2 > the operator (Op^Jcr])^ 1 ((Op h [cr])*) £2 belongs to the trace class and 

(2.6) l\m(2nh) d - 1 Tr{(Op h [a])^{(Op h [a]yY 2 )= [ [ a (u , O^^W^du . 

Sketch of proof. First suppose that m > d — 1 and £\ = 1, £2 = 0. Then Op/Jcr] is trace 
class and the asymptotics (12.61) follow by substituting the representation ( 12. 5ft into ( 12. 4p . 
In the general case, using local coordinates on the sphere, iterating a standard composition 
formula for \l/DO in L 2 (IR d_1 ), and using the formula for the symbol of the adjoint operator, 
we obtain that for any N > 0, 

(Op>r((Op h [a])f 2 = Pfe [x] +R N (h), 

where: 

(i) the symbol x G S((£)~ m£l ~ mt2 ) can be represented as 

x = x + hx\ 

with x , Xi G S((^)~ m£l ~ m£2 ), x is independent of h and 

(ii) -Rtv(^) is an operator with the integral kernel in C^S^ -1 x § d_1 ) and the C^-norm of 
R N (h) is 0(h N ) as /1 0. 

This reduces the problem to the case £\ = 1, £2 = 0. □ 

Remark. In the situation we are interested in, the representation (12. 5p arises as a corollary 
of the asymptotic expansion 

00 

j=0 

but we are only interested in the first term of this expansion. 

In our construction, the \l/DO will be defined in terms of their amplitudes rather than 
their symbols. Thus, we need a statement which is standard in the \l/DO theory (see e.g. 

TO- ' 

Proposition 2.2. Let m > 0, and let b = b(u, u', £, h) be a smooth function of the variables 
(w,£) G T*S d_1 ; u' G § d_1 and h G (0, 1). Assume that b satisfies the estimates 

(2.7) \d£dZeZ,b{u>,u/,t,h)\ < CVy(£)- m -H 

for all multi-indices a, (3, 7 uniformly in h G (0, 1) . Then for any N > 0, the operator with 
the integral kernel 

(2.8) (2nh)- d+1 [ e-'^-^ )/h &(w,c/,£,/i)d£ 

can be represented as Op h [<r] + R]y(h), where the following conditions are met: 
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(i) The symbol a can be written as a = cr + ha\ with cr , o\ G <S((£) _m ) and 

ao(w,€,h) = b(w,w,£,h). 

(ii) The operator i?jv(^) has the integral kernel in C^S^ -1 x with C N -norm satisfying 
0{h N ) ash^O. 

2.3. Approximate solutions to the Schrodinger equation. Here we recall the con- 
struction of approximate solutions to the Schrodinger equation Hu = k 2 u from [17]. The 
solutions u are sought as functions 

u = u(x,p), x G M. d , p<ER d , \p\ = k. 

We denote p = plpl" 1 G S> d_1 . Writing 

u(x,p) = e ie ( x,p ^v(x,p), 

we obtain the eikonal equation 

|V0| 2 - 2(A, V8> + (V(x) + \A{x)\ 2 ) = k 2 
for the phase function and the transport equation 

-2z(V6, Vv) +2i(A, Vv) - Av + {-iAQ + idivA)v = 

for the amplitude function v. The approximate solution to the eikonal equation is constructed 

as 

(2.9) @±(x,p) = (x,p) +(p±(x,p), 

POO 

(2.10) (j) ± ( x ,p) = T (A(x±tp),p)dt. 

Jo 

Next, for a given iV G N, the approximate solution to the transport equation is con- 
structed as the asymptotic series 

N 

(2.11) v^\x,p) = ^(2ik)- n v^(x,p), 

n=0 

where Vq(x,p) = 1 and the coefficients Vn are determined via an explicit iterative proce- 
dure: 

„,(±) 



/■oo 

\(x,p) = =F / fi ±) (x±tp,p)dt, 
Jo 



v (±) 



where 

= 2*(A - V0±, V^) - Av^ + (| V0±| 2 - 2(A, V0±) + V! - zA0 ± ; 
and 

Vi = |A| 2 + idivA 
The functions </>± and solve the eikonal and transport equations up to error terms that 



can be explicitly controlled; we do not need the precise statement here, see [TyJ for the 
details. Here we need <f>± and v± simply as "building blocks" for the approximation to the 
scattering amplitude, which is given in the next subsection. 

When considering the functions (f>± and Vn , we will always exclude a conical neighbour- 
hood of the direction x = —p (for the sign "+") or x = p (for the sign "— "). Outside these 
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neighbourhoods, the functions <p± and v^t\ n > 1, decay at infinity in the x- variable. More 
precisely, the following statement is proven in [17] : 

Proposition 2.3. [17J Let assumption (11.11) hold and let x G M. d , uj G S d_1 be such that 
±(x,uj) > — 1 + e for some e > 0. Then the functions (j)± and Vn , n > 1, satisfy the 
estimates 



(2.12) \d a XM*^)\<C aP (x) l -p-\«\ 



(2.13) \d:d^\x,u)\<C a ^(x)- n ~\ a \, 
for all multi-indices a, (3. 

We will write 

(2.14) u { ±\x,p) = u ± (x,p) = e ie±{x > p) v ( ±\x,p). 

2.4. Approximation to the scattering amplitude. Here we recall the approximation 
to the scattering amplitude obtained in [T7J. It is known (see pQ) that off the diagonal 
uj = uj' , the integral kernel s(uj,uj'; k) of the scattering matrix S(k) is a C^-smooth function 
of co, uj' G and it tends to zero faster than any power of k~ x as k — \ oo. Thus, it suffices 
to describe the structure of s(u>, u'; k) in a neighbourhood of the diagonal uj = uj'. Fix some 
5 G (0, 1); for an arbitrary point ujq G S d_1 , let fl(uJo) C S d_1 be the conical neighbourhood 
of ujo given by 

(2.15) n(u ) = {uj G S d_1 : (to, coo) > 5}. 
Let u± be as in (I2.14p . We set 

duj u = (Vm,w ), 

where the gradient of u = u(x,p) is taken in the x variable. For uj,uj' G Q(ujq), define 

(2.16) s { N \uj,uj';k) = -ink d - 2 (2n)- d x 

dx— 



x 



u^\x, kuj)(d LUO u ( L S) )(x, kuj') — (d Uo u+)(x, kuj)u^\x, kuj') 



2i / (A(x),uj }uf\x, ku)u^\x, kuo')dx . 



The integrals in (I2.16P do not converge absolutely and should be understood as oscillatory 
integrals. In other words, (12. 16[) should be understood as a distribution on Q(ujo) x Q(uj ). 

Proposition 2.4. [TTJ For any q G M there exists N = N(q) G N such that for any ujq G 
the kernel 

s^- 1 (uj, uj'; k) = s(uj, uj'; k) — (uj, uj'; k) 



belongs to the class C q (fl(u ) x Q(uj )), and its C q -norm is 0{k q ) as k — >■ oo. 
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2.5. The scattering matrix as a $DO on the sphere. Below we represent the scattering 
matrix S(k) as a semiclassical \l/DO on the sphere. The semiclassical parameter is h = fe -1 . 
The idea to represent the scattering matrix as a \l/DO on the sphere goes back to |4J. The 
statements almost identical to Lemma 12.51 can be found in [T7J Propositions 6.1 and 6.4] 
and [121 Section 8.4], but for technical reasons these statements are not sufficient for our 
purposes. A related work (but written from a very different viewpoint) is [3J, where the 
scattering matrix is represented as a Fourier integral operator corresponding to the classical 
scattering relation. 

Lemma 2.5. Let assumptions fll.l[) hold, and let m = min{l,p — 1}. For any q G N, the 

scattering matrix can be written as 

(2.17) S(k) = I + Pk - 1 [a}+R g (k), 
where: 

(i) The symbol a can be represented as 

(2.18) <7 = <7 + ATVi, 

with (To, 0i G iS((£)~ m ) and 

(2.19) < 7 (w,0 = exp(iM(a;,0)-l; 

(ii) the operator R q {k) has an integral kernel in the class C^S ' -1 x S d_1 ) and its C q -norm 
is 0{k~ q ) as k — )■ oo. 

Proof. 1) Let ipi,ip2 £ C 00 ^ -1 ) be functions with disjoint supports. Then ^iS(k)ip2 has a 
C^-smooth integral kernel which decays faster than any power of fe _1 as k — >■ oo. The same 
comment applies to ipi Op fc -i [a}ip 2 with a G S({x)~ m ). This shows that using a sufficiently 
fine partition of unity on the sphere, one easily reduces the problem to approximating the 
integral kernel of S(k) locally in any conical neighbourhood Q(uq), see f)2.15p . Thus, we can 
use Proposition 12.41 

2) Let us rearrange the integrand in f)2.16p . Denote 

(2.20) w±{x,p) = e i<t>±ix '®v { ±\x,p), 

(2.21) w±(x,p) = ke i<t,±{x 'P\v ( ±\x,p) - 1), 
so that from (12. 14ft . 

u±(x, kuj) = e ik{x ^w ± {x, few) = e ik M^±M + k~ 1 w±(x, few)), 

(d Uo u±)(x, few) = e^ x ^[ik(u ,u})e^ x ^ + z(w , w)w±(z, few) + (^w±)(x, few)]. 
Now some elementary algebra shows that formula (I2.16P can be rewritten as 

(2.22) s { "\u,u';k) = (J^j J '^a(w, w', x)dx, 
where 

(2.23) a(w, w', x) = -(w , w + w') exp(i0_(x, w') — i<fi+(x, w)) + fe _1 ai(w, w', x), 
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(2.24) 2 ai (u,u',x) = (u, u ) (e-^^uLfo ku') + w+(x, ku)w-(x, few')) 
+ (u/, w ) (e^-^'^w+ix, ku) + w + (x,ku)w„ (x, ku')) 



+i(d Uo w + )(x, ku)ui-(x, ku')—iw + (x, ku)(d OJo W-)(x, ku')— 2(A(x), uq)w + (x, ku)w-(x, ku'). 

Note that the choice a(u, u', x) = ^(u + u', u ) in (12.221) yields a 5-function on the sphere. 
Thus, we can write 

(2.25) s { N \u,u';k) - 5(u - u') = (— \ [ e^"-"' ' x \a + kr l ai ){u, u', x)dx, 
where 

(2.26) ao( w , x ) — 2 ( u o> u + u ') (exp(i0_(x, u') — i(f)+(x, u)) — l). 

3) Let us change variables in the integral (12.251) . Instead of integrating over x G A Wo , we 
shall integrate over (6 A u , where 

(2.27) x = £- ^ > + o/). 

(tu + w ,a;o) 

Recall that u,u' G fi(wo), and so the denominator in (12.271) does not vanish. An inspection 
shows that 

(x, u — u') = (£, u — u'). 

Thus, we obtain 

(2.28) s^\u,u';k)-5(u-u')^ (J^J *jf e- ik ^'^b(u,u',0^, 
where b = bo + A; _1 6i with 

(2.29) bj(u,u',£) = J(u,u')aj(u,u' ,#(£)) 

and J(u,u') denotes the Jacobian of the linear map (I2.27j) considered as a map from A w to 
A Wo . It is easy to see that J(u,u') is a smooth function of w,u/ G f2(a;o) and 

(2.30) J(w,w) = (t^Wo) -1 . 

4) The r.h.s. of (I2.28P is a semiclassical \l/DO with the amplitude 6 and a semiclassical 
parameter /i = A; -1 , see (12. 8p . In order to complete the proof, by Proposition 12.21 if suffices 
to check the estimates 

(2.31) < cW£r mHQ| , 

for j — 0, 1 and all multi-indices a, /3, 7 uniformly over k > 1, and to check the identity 

(2.32) 6 (w, w, = exp(zM(w, 0) - 1. 

Let us first check (12.321) . Recalling the definition (I1.3P of M and the definition (I2.10p of (fi±, 
we get 

From this and fl2T26|) . fl2T29"l) and ( 12730]) . we obtain 

6o(^,o;,0 = J(w,w)ao(w,t«;,a;(0) = exp(iM(o;, x(£))) - 1, 
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where is the linear map (12.27)) . Next, by the definition of the map x(£), for uj = oj' it 
takes the form x(£) = £ + ecu, and by the definition of the function M we have 

M{u,£ + cu)=M{u,0- 

Thus, we obtain (12. 32 p . 

5) It remains to check that the estimates ( 12. 31 )) are satisfied. This is essentially a consequence 
of Proposition |2T3| let us check this. Recalling that m = min{l, p— 1} and using the estimates 
( 12T2|) . (12TT3|) . we obtain 

\d^Mx,u)\<C a ,(x)- m -^, 

\d:d^w ± (x,ku)\<C aP (x)- m -^, 

\d%di{d Uo w ± ){x,ku)\ < C aP (x)- m -\ a \ 

\d^w ± (x,koj)\<C aP (x)-^, 

where all the estimates are uniform in k > 1. It follows that Oq and a 1; defined by ( I2.26j) . 
( I2.24p respectively, satisfy 

(2.33) |^d>>,u/,*)| < C aP (x)- m -\ a \ 

uniformly in A; > 1. Now from (12.271) . (12.29)) and (12.331) by an elementary calculation we 
obtain ACT]) . □ 

3. Proof of Theorem 11.11 
3.1. The case of a monomial (p. 

Lemma 3.1. Assume the hypothesis of Theorem \l.l\ Then for any integers l\ > 0, £2 > 
such that i\ J ri2>id—V)jm,m = min{l, p — 1}, the asymptotic formula fll . 16j) holds true. 



Proof. By Lemma [2.51 we have 

(3.1) (s(k) - iy>(s(ky - if* = (Op^*] + R q (k))^ ((Op^M)* + R q (ky) e \ 

where a, R q (k) are as described in Lemma [2.51 Expanding the brackets in (13. ip . we obtain 

(3.2) (S(k) - I)^(S(ky - It = (Ov k ^]) h ((Ov k -^\rf + Q q (k), 

where Q q (k) is the sum of the products of operators, to be estimated below. For the first 
term in the r.h.s in ( 13. 2p . by Proposition 12.11 we have 

Let us check that by a suitable choice of q we can ensure that the error term Tr Q q (k) remains 
bounded as k — > 00; this will certainly yield the desired asymptotics ( 11.161) . By choosing q 
sufficiently large, we can make sure that the estimate (||-||i is the trace norm) 

(3.3) 11^)111 = 0(1), k ^00, 
holds true. Next, using the estimates 

|Tr(AB)| < HAHH^Hi and ||0|| < ||C||i 
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and recalling that Q q {k) arose as a remainder term in the expansion of the brackets in the 
l.h.s. of (13. ip . we obtain 

(3.4) \Tr(Q q (k))\ < C(^i,4) max {\\R q (k)\\{,\\Op k -i[a}\\ il+i2 ' j }. 

By (12. 3p . we have 

IIOp^-xHH = 0(1), k^oo. 

Combining the last inequality with (13. 3p . we obtain that Ti(Q q (k)) is bounded as k — > oo, 
as required. □ 

3.2. Application of the Weierstrass approximation theorem. 

Lemma 3.2. Let £$ be an even natural number and let v, v k (k > I) be a-finite measures 
on T \ {1} such that 



(3.5) 



/ \z - l\ £ °dv{z) < oo, f\z- l\ e °du k {z) < oo, 
Jt Jt 

for all k. Suppose that for all integers £\ > 0, £2 > such that £\ + £2 > £0, the relation 
(3.6) lim f(z- l) ei (z- lf 2 du k {z) = f (z - l) ll {z- \) h dv{z) 

k^oo J j J j 

holds true. Then for any ip G C(T) such that f(z)\z — 1| is continuous on T, the relation 



(3.7) lim / ip(z)dv k (z) = j (f(z)di i (z). 
holds true. 

Proof. Applying the Weierstrass approximation theorem to the function f(z)\z — 1| °, for 
any e > we obtain a polynomial ipo{ z ) m z, z such that 

\<p(z)\z-l\- e °-<p (z)\<e, feel 

Let us define (f±(z) = (Re (po(z) ± s)\z — then it follows from the above that 

(3.8) ip_{z) < Re<p{z) < ip + (z) Vz G T, 

(3.9) <p + (z) -<p-(z) = 2e\z- l\ e °. 

By the construction of <p±, it can be represented as a polynomial 'm w — z — 1, w = z — 1 
involving only products w ll w l2 with £\ + £ 2 > £q- Thus, by (13. 8ft . (13. 9 p we can write 



(3.10) / y>_(z)dz/(z) < / Re(p(z)dv(z) < / y? + (z)di/(z), 
Jt Jt Jt 

(3.11) / if-{z)du h (z) < / Re^(z)dz/ fe (2) < / (p + (z)du k (z), 



(3.12) / p+(*)dz/(;s) - / (p_(z)djs(z) =2e \z - l\ e °du(z), 

Jt Jt Jt 
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where all integrals are absolutely convergent by (13.51) . Now we can use (13. 6p to pass to the 
limit in (13. lip . Using (I3.10p . (I3.12p and denoting by C the value of the integral in the r.h.s. 
of (I3.12p . we obtain 

limsup / Re (p(z)dv k (z) < / (p + (z)du(z) < / Re(p(z)du(z) + 2eC, 

liminf / Re (p(z)dvk(z) > / (f^(z)du(z) > / Re (p(z)du(z) — 2eC. 
Since e > may be taken arbitrary small, this yields 

lim / Re Lp(z)dh>k(z) = / Re (p(z)dv(z). 



Since the same argument can be applied to the imaginary part of ip, we obtain the required 
statement. □ 

3.3. Proof of Theorem II. li Recall that the measures /i^ and /i are defined in ( II. 5p and 
( 11. 6p . By (II. 7p . ( II .8p . the conclusion of Lemma I3TT1 can be written as 

lim k- d+1 f(z- l) h (z-lY*d^ k (z) = f(z - l) h (z-lY*dfi(z). 

Now it remains to apply Lemma [3.21 with v k = k~ d+1 ^ and v = fi. 
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